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Abstract. Group classification of a class of Benjamin-Bona-Mahony (BBM) equations with 
time dependent coefficients is carried out. Two equivalent lists of equations possessing Lie 
symmetry extensions are presented: up to point equivalence within the class of BBM equations 
and without the simplification by equivalence transformations. It is shown that the complete 
results can be achieved using either the gauging of arbitrary elements of the class by the 
equivalence transformations or the method of mapping between classes. As by-product of 
the second approach the complete group classification of a class of variable-coefficient BBM 
equations with forcing term is derived. 


1. Introduction 

The third-order nonlinear partial differential equation 

Ut T Ux T UUx '^xxt — 0? 

named these days the Benjamin-Bona-Mahony (BBM) equation, appeared in [T] as an 
alternative to the Korteweg-de Vries equation, ut + Ux + uux + Uxxx = 0, model for the 
unidirectional propagation of moderately long waves with small but finite amplitude in systems 
that manifest nonlinear and dispersive effects. Numerical studies showed that the BBM equation 
admits soliton solutions whose interaction is inelastic though close to elastic m- It was proved 
in [3] for the equivalent form ut = uUx + Uxxt of the BBM equation that it has no conserved 
quantity in addition to those found by Benjamin, Bona and Mahony: u (mass), {v? + u^)/2 
(energy), and m^/ 3 (momentum). Lie symmetries and the corresponding reductions of the BBM 
equation in the above equivalent form were obtained in [5] (these results were also presented 
in [6l pp 194-196]). It was found that the maximal Lie symmetry algebra of this equation is 
a three-dimensional Lie algebra of the type A 2.1 0 Ai spanned by the vector fields dt, tdt — udu 
and dx- 

There are several recent works (see [7j and references therein) devoted to the study of Lie 
symmetries of variable-coefficient BBM equations of the general form 


Ut + f{t)ux + g{t)uux + h{t)Uxxt = 0, 


( 1 ) 


where /, g and h are arbitrary smooth functions of the variable t with gh ^ 0. However none 
of these works contains exhaustive and completely correct results. We aim to fill up this gap by 
presenting the exhaustive group classihcation of equations from class dU) and classifying local 
conservation laws of these equations. 

In order to study transformational properties of class in Section 2 we describe the entire 
set of its admissible transformations (i.e., the equivalence groupoid of this class) by proving that 
the class is normalized in the usual sense. The assertions concerning reducibility of variable- 
coefficient BBM equations to their constant-coefficient counterparts are presented therein. The 
complete classification of Lie symmetries is obtained in Section 3 under gauging the arbitrary 
element 5 to 1 by a parameterized family of point transformations related to the associated 
equivalence group. The method of mapping between classes is utilized in Section 4 to verify 
the classification results. As by-product of this approach the complete group classification of 
a related class of variable-coefficient BBM equations with forcing term, ut+uux + h{t)uxxt = 
is derived. In Section 5 we construct local conservation laws of equations from the class ©• 

2. Admissible transformations 

Two differential equations which are connected by a point transformation are called similar [5j. 
Similar differential equations have similar sets of solutions, symmetries, conservation laws 
and other related objects. Therefore, it is instructive to study point transformations linking 
equations from a given class of differential equations. Such transformations are called form¬ 
preserving |16] or allowed [18]. The rigorous reformulation of this notion leads to the notions 
of admissible transformations laiinj and the equivalence groupoid of the class mm- An 
admissible transformation is a triple that consists of two similar equations from the class and 
a point transformation between them. The set of admissible transformations of the class possesses 
the groupoid structure with respect to the natural operation of the composition of admissible 
transformations and hence it is called the equivalence groupoid of the class. 

A joint transformation of the unknown functions and the arbitrary elements is called 
an equivalence transformation of the class if it satisfies the following properties: (f) The 
transformation becomes point with respect to the independent and the dependent variables 
when the arbitrary elements are fixed, [ii) It is consistent with the contact structure of the 
underlying jet space. [Hi) It maps any equation from the class to an equation from the same class. 
Depending on additional restrictions, one have different kinds of equivalence transformations 
(usual, generalized, extended and generalized extended ones). Equivalence transformations of 
the same kind constitute a group, which is called the equivalence group of the class (of the 
respective kind), see [TU1[T3| and references therein. If the equivalence groupoid of the class is 
generated by the equivalence group of this class, then this class is called normalized in the same 
sense as the kind of the equivalence group involved mm- 

One-parameter subgroups of the equivalence group can be found by the Lie inhnitesimal 
method, whereas the direct method [SHIS] allows one to construct the entire equivalence 
group including discrete equivalence transformations and this technique is preferable. The direct 
method can also be used for finding the equivalence groupoid and checking whether the class 
is normalized [TOHSHIT]. A very useful feature of normalized classes is that the equivalence 
groups for their subclasses, singled out by setting additional constraints on arbitrary elements, 
are subgroups of the equivalence group of the entire class. We use this property to derive the 
equivalence group of the class ([T]) from the equivalence group of a normalized superclass. 

Consider the class of (l-l-l)-dimensional partial differential equations of order p + q > 1 with 
hxed p,q ^ No, 

Upq H (f, 1)) ) (2) 

where Upq = d^^’^u/dt^dx'^ and _i) denotes the set of all derivatives of u up to order p+q— 1 


including u as the zero-order derivative. The function H is an arbitrary smooth function of its 
arguments. The following assertion was proven in m- 

Theorem 1. A point transformation t = T{t,x,u), x = X{t,x,u), u = U{t,x,u), where 
\d{T, X,U)/d{t,x,u)\ 7 ^ 0, maps a fixed equation from class ([2]) to an equation from the same 
class if and only if it satisfies the following conditions: 

(i) T = T{t), X = X{x) for pq^O, p^ q; 

(ii) T = T{t), X = X(x) or T = T{x), X = X{t) for p = q^ 0; 

{in) T = T{t) for p = 0, q 0; 

(iv) X = X{x) for p 0, q = 0. 

Corollary 1. The usual equivalence group of class ([2]) with any fixed p and q consists of 

the transformations whose projections to the space of {t, x, u) coincides with that described in 

Theorem 1. Moreover, this class is normalized in the usual sense. 

Therefore, Theorem 1 provides the complete description of admissible (form-preserving) 
transformations in class m- As this class is normalized, we can use transformation constraints 
from Theorem 1 for describing admissible transformations in subclasses of this class, in particular 
in class ©• 

For class ([T]) we have p = 1 and q = 2, which corresponds to Case (i) of Theorem 1. Now we 
directly seek for all point transformations of the general form 

i=T{t), X = X{x), u = U{t,x,u) (3) 

with TtXxUu 7 ^ 0 that map a fixed equation of the form ([T|) to an equation from the same class, 

% + fit)ux + g{i)uux -t- h{i)u^^i = 0. (4) 

The components of the prolongation of transformation ([3]) to the involved partial derivatives are 

hi — 7=- {Ut T bJuUt) , Ux — „ {Ux Uu'^x) j 

hxxi - rp (fhtx T UxuUf T Utu'^X T bJ^tu'^t'^x T UyUtxf) ~\~ 

T rr\ ta2 T UxxuUt T 2Ufxu'^X T 2Uxuu'^t'^X T 2Uxu'^tX T bjfuu'^x T bJuuu'^t'l^x T 

It^X 

T 2UuuUxUtx T UtuUxx T UuuUfUxx T bJuUxxt) ■ 

Using these expressions, we rewrite equation (jH) in the old variables, which are without tildes. 
Then we substitute Uxxt = —{ut + f{t)ux + g{t)uux)/h{t) into the obtained equation in order 
to conhne it to the manifold dehned by ([T]) in the third-order jet space with the independent 
variables (t, x) and the dependent variable u. Splitting the resulting identity with respect to 
the derivatives ut, Ux, utx and Uxx leads to the determining equations on the transformation 
components T, X, and U, 


Uuu-Utu-b), XxxUu - 2XxUxu, XxUu - h{XxxUxu XxUxxu) + ^XxUu, 


Ut + ^{f + gU)Ux = h^Utx - h^Uxxt, f + ~gU = 


The equations Uuu = Utu = 0 imply the representation U = U^{x)u + U^{t,x), where and 
are smooth functions of their arguments with 7 ^ 0. We substitute this representation for U 






into the remaining determining equations and further split with respect to u. As a result, we 

get the equations = Ul = = 0, / = - U°g), g = and h = X^h. 

We solve the equations for X, and U^, and substitute the obtained expressions into the 
equations representing transformation components for arbitrary elements, which leads to the 
following assertion. 

Theorem 2. The usual equivalence group Gi of class © is comprised of the transformations 
t = T{t), x = Six + 62 , u = S3U + 64, f = - ^4g), g = 7^9, h = 5lh, 

dt03 dtOs 

where 5j, j = 1, 2, 3,4, are arbitrary constants with ^ 0 and T = T{t) is an arbitrary smooth 
function with Tt ^ 0. Class © is normalized in the usual sense. 

Thus, each point transformation between equations from class ([T]) is induced by an element 
of the group . In order to find which variable-coefficient equations of the form ([T]) admit 
constant-coefficient counterparts, we assume that the transformed arbitrary elements /, g and h 
are constants in equivalence transformations. This results in the following assertion: 

Proposition 1 . A variable-coefficient equation from class ([T]) is reduced to a constant-coefficient 
equation from the same class by a point transformation if and only if the corresponding 
coefficients f, g and h satisfy the conditions 



i.e., h is a constant and f is proportional to g. 

Equivalence transformations allow us to simplify the consideration by reducing the number 
of arbitrary elements. For example, we can set the gauge g = 1 using the family of point 
transformations 



parameterized by the arbitrary element g and related to equivalence transformations from the 
group Gf. Then the other arbitrary elements are f{f) = f{t)/g{f), and h{t) = h{t). Here and 
below an integral with respect to t should be interpreted as a fixed antiderivative. 

Therefore, without loss of generality we can restrict ourselves by the study of the class 

ut -\- f{t)ux -C uux -h h{t)uxxt = 0, (6) 

as all results on symmetries and exact solutions for equations of the form ([ 1 ]) can be recovered 
from the corresponding results obtained for equations of the form ([ 6 ]) using the above family of 
point transformations. 

Since class © is normalized, the equivalence group of its subclass ([ 6 ]) can be easily found as 
the subgroup of the group Gf whose elements preserve the gauge <7 = 1. 

Corollary 2. Class ([ 6 ]) is normalized in the usual sense. Its usual equivalence group Gf is 
constituted by the transformations 

i=^t + 6o, x = 5ix + 52, u = 63u + 5i, f = dsf — 84, h = dfh, 

03 

where j = 0 ,...,4, are arbitrary constants with ( 5 i (53 7 ^ 0. 

There are no truly variable-coefficient equations in class ([ 6 ]) that are reduced by point 
transformations to constant-coefficient equations from the same class. 

In the next section we use the equivalence transformations in the course of group classification 
of classes © and dH). 






3. Lie symmetries 

We carry out the group classification of class ([ 6 ]) within the framework of the classical Lie 
approach mm- We look for vector fields of the form 

r = T{t, X, u)dt + X, u)dx + r]{t, x, u)du 

which generate one-parameter groups of point-symmetry transformations of an equation from 
class Q and hence jointly constitute the maximal Lie invariance algebra of this equation. 
Any such vector held, L, satishes the inhnitesimal invariance criterion, i.e., the action of the 
third prolongation, of the vector held F on equation ([ 6 ]) results in the conditions being an 
identity for all solutions of this equation. That is, the criterion is read as 

r(3) [ut + f {t)u^ + UUx -b h{t)Uxxt] \u^^t=-{ut+f{t)u^+uuA/Ki) = (7) 

Theorem 2 allows us to simplify the computation since it implies that r = T(t), ^ = ^(x), 
7] = r]{u), and = Vuu = 0, Then the left hand side of the condition ([7]) is a polynomial in the 
variables x, u, ut and Ux with coefficients depending on t. Hence these coefficients are equal to 
zero, which gives the determining equations on remaining arbitrarinesses in the components r, 
^ and r] of the vector held T. Solving the simplest of these equations we get the general form of 
the components of a Lie symmetry operator of an equation from class (j 6 |l. 

T = (^C2 - Cl) t-b Co, C = C2X-bC3, 7 / = (ic2 "b Cl) U-b C4 

with Cj, f = 0,..., 4, being arbitrary constants. The residuary determining equations involving 
arbitrary elements have the form 

{{^C2 - Ci)t + Co)ht = 2 c 2 h, ((^C2 - Cl) t-b Co)/i = (^C2 -b Cl)/ - C4 ( 8 ) 

and are called classifying equations. Further analysis is carried out using the method of furcate 
splitting introduced in [22]. For each vector held F from the algebra equations (| 8 |) gives 

a system on / and h of the general form 

{at + j3)ht = -fh, {at + P)ft = {^-y - a)f + u, (9) 

where a, /3, 7 and zz are constants. Varying vector helds from we obtain a set of systems of 
the form ([S]). The number k of such systems with linearly independent tuples of the coefficients 
a, /3, 7 and v is not greater than two; otherwise the set of these systems would be inconsistent 
in total. Therefore, possible values of k are /c = 0, /c = 1 and k = 2. 

If A; = 0, then equations dSj) are identities with respect to / and h and hence cq = ci = C 2 = 
C 4 = 0. Thus, we obtain that the maximal Lie invariance algebras of equations from class dbj) for 
general values of / and h is the one-dimensional algebra {dx), which gives Case 0 of Table 1. We 
get the same result when varying the arbitrary elements / and h an splitting equations dHj) with 
respect to them. In other words, the kernel of the maximal Lie invariance algebras of equations 
from class ([ 6 |) (i.e., the common part of the maximal Lie invariance algebras of equations from 
this class, see | 8 | for the dehnition) coincides with the maximal Lie invariance algebra {dx) of 
the general case. 

If /c = 1, equations ([9|) are not identities but equations on / and h. In order to get Gif- 
inequivalent solutions of these equations, we study how transformations from Gf act on the 
coefficients a, /?, 7 and v of system (l9|). Combined with the multiplication by a nonzero 
constant y, each transformation from the equivalence group Gf is extended to the coefficient 
quadruple of the system (l9|) as 

a = ax, ^ 7 = 7 X, k = {63U + 64 {^J - a)) x- 

Analysis of the induced transformations for a, /?, 7 and n results in the following assertion. 


Lemma 1. Up to G 2 -equivalence the parameter quadruple (a,/?, 7 , 1 ^) can be assumed to belong 
to the set 

{(1,0,;9,0), ( 1 , 0 , 2 , 1 ), ( 0 , 1 , 1 , 0 ), ( 0 , 1 , 0 ,!)}, 
where p is an arbitrary constant. 

The group classification for the case A: = 1 is obtained as follows. Seqnently solving system ([9]) 
for each of the four inequivalent valnes of (a, /3, 7 , n) presented in Lemma [U we get the 
corresponding cases for the arbitrary elements {h,f), 

[et^,Xt^), (et^,lnt + A), (ee^Aea), (e,t + A), 

where e and A are constant parameters arising in the course of integration, and e 7 ^ 0. Then we 
check whether the remaining parameters can be additionally gauged by transformations from 
the group It appears that s can be scaled to 1 if it is positive and to —1 if it is negative. 
In the second and the fourth cases the constant A can be set eqnal to zero by a translation of u. 
In the first and the third cases A cannot be gauged. We substitute the obtained forms of h and / 
into equations ([5|) to get constraints on constants c*, i = 0,..., 4. This gives the respective forms 
of Lie symmetry operators. The results are collected in Cases 1-4 of Table 10 

If /c = 2, then up to G^-equivalence and linear combining there exists a single pair of the 
coefficient tuples, namely, ( 1 , 0 , 0 , 1 ) and ( 0 , 1 , 0 , 0 ), for which the corresponding systems of 
the form Q are consistent to each other. Then / and h are constants, / = 0 mod and 
h = e mod G^, where e = ±1, which gives Case 5 of Table 1. 

We sum up the above consideration. 

Theorem 3. The kernel of the maximal Lie invariance algebras of equations from class (l 6 |) is 
the one-dimensional algebra (dx). All possible Gf-inequivalent cases of Lie symmetry extensions 
are exhausted by Cases 1-5 of Table 1. 

Corollary 3. The group elassification for class © up to GI-equivalence results in the list 
presented in Table 1, where the arbitrary element g is assumed to be equal 1. 

In order to get the classification list for class ©, where forms of arbitrary elements are 
not simplified by equivalence transformations, we apply transformation © combined with 
transformations from the equivalence group Gf to equations of the form Q with / and h 
presented in Table 1. Basis elements of the corresponding maximal Lie invariance algebras are 
pushed forward by the same transformations. Then we re-denote the constants and collect the 
obtained results in Table 2. The detailed procedure of the equivalence based approach for deriving 
most general forms of arbitrary elements and basis elements of the corresponding maximal Lie 
invariance algebras can be found in m- 

4 . Mapping between classes and group classification of the related class 

An alternative way for group classification of class ([T]) is the method of mapping between classes, 
which was suggested in [23]. This method has been successfully applied to several classes of 
nonlinear partial differential equations (see, e.g., |2Ujl. 

Class ([I]) can be mapped to a similar class of third-order partial differential equations of the 
form 

Ut uUx -\~ hi^t^Uxxt — ^(A), h 7^ 0. (lb) 

^ Arguments of logarithm and bases in powers with noninteger exponents containing p should be carefully treated 
throughout the paper. In fact, one should take their absolute values if they can be negative. Usually the positivity 
of t well agrees with the physical meaning of t as time variable. So, omitting the notation of absolute value is not 
too essential in Table 1. At the same time, the expression eT + k from Table 2 may be negative and hence we 
implicitly assume replacing it by its absolute value. 


Table 1. The group classification of class ([6]) up to -equivalence. 


no. 

h{t) 

fit) 



Basis of 

0 

V 

V 

dx 



1 

etP 

p —2 

Xt^~ 

dx, 

2td 

t + pxdx + ip - 2)udu 

2 


Int 

dx, 

tdt 

1 

+ 

3 

ee^ 

Ae^* 

dx, 

2dt 

-1- xdx + udu 

4 

e 

t 

dx. 

dt- 

- du 

5 

e 

0 

dx. 

dt. 

tdt - udu 


Here p and A are arbitrary constants, e = ±1 mod G^- In Case 1 (p, A) / (0,0). 


Table 2. The group classification of class m using no equivalence. 


no. 

h{t) 

fit) 


Basis of 

0 

V 

V 

dx 


1 

piieT + k)p 

P-2 

/i25(eT + K) ^ + p3g 

dx. 

2 

-{sT + K)dt + epxdx -b e{p - 2){u + P 3 )du 

9 

2 

piieT + 

P2fl'ln(eT + k)+ P39 

dx, 

-{eT + n)dt + exdx - ep 2 du 

9 

3 

pi exp((TT) 

P25exp(lcjr) + p3g 

dx. 

2 

-dt + axdx + (t{u + 93)du 

9 

4 


P 29 T P39 

dx, 

~dt f^2du 

9 

5 


939 

dx, 

1 T 

-dt, —dt -iu + p3)du 

9 9 


Here g is an arbitrary nonvanishing smooth function, T = J’p(f) df; e = ±1; pi, /i 2 , Ps, and p 
are arbitrary constants satisfying the following constraints: piX ^ 0; in Case 1 pp 2 7^ 0; and in 
Case 4 p 2 7 ^ 0. 


The map is realized by the family of point transformations 


t = 



X = X, 


u = u + 


m 

9{t) ’ 


( 11 ) 


parameterized by two arbitrary elements of class m- The arbitrary elements in the imaged 
equations take values (tildes in (fTn|) are omitted) 




Following the method of mapping between classes, we first classify Lie symmetries of the imaged 
class (jlOp and then use the family of point transformations (HID to extend the result to the initial 
class (dl). 

In order to efficiently solve the group classification problem for class m, we look 
for admissible transformations in this class using the direct method. It appears that such 
transformations are exhausted by transformations from the usual equivalence group of this class. 











Theorem 4. The usual equivalence group of class cni) consists of the transformations 


~ 

t — —t + 6 o, 

03 


X = 5ix + 82 , u = 63 U, h = 5\h, 


l=^Il 


where 5j, j = 0 ,1,2,3, are arbitrary constants with JiiJa 7 ^ 0 . Class (fT 0 |) is normalized in the 
usual sense. 

Using the classical Lie infinitesimal method in the same way as in the previous section, we 
get the complete group classification of equations from class (llOp . The results are summarized 
in the following assertion. 

Theorem 5. The kernel of the maximal Lie invariance algebras of equations from class (unD is 
the one-dimensional algebra (dx). All possible -inequivalent cases of extension of the maximal 
Lie invariance algebras are exhausted by Cases 1-4 of Table\^ 


Table 3. The group classification of class (fT 0 |) up to G^-equivalence. 


no. 

h{t) 

l{t) 


Basis of 

0 

V 

V 

dx 


1 

£tP 

Xt^ 

dx, 

2 tdt + pxdx + (p - 2 )udu 

2 

ee^ 

Aei* 

dx. 

2 dt + xdx + udu 

3 

e 

1 

dx. 

dt 

4 

£ 

0 

dx. 

dt, tdt - udu 


Here p and A are arbitrary constants, e = ±1 mod Gf, in Case 1 (p. A) 7 ^ (Q, 0). 


Remark 1. The most general forms of the functions h and I that correspond to equations from 
class m with Lie symmetry extensions are 

1 . h = Xi{et + k)p, I = X 2 {st + k)^: = (^Qx, 2{et + K)dt + spxdx + e{p — 2)udu); 

2 . h = Xie^^, I = A 2 e 2 '^*: 20 ^ + axdx + crudu)] 

3. /i = Ai,/ = A2: A^^ = {dx,dt); 

A. h = Xi, I = 0 : A'^^^ = {dx, dt, tdt - udy). 

Here Ai, A 2 , £, k and p are arbitrary constants with Xiae 7 ^ 0. Additionally, in Case 1 
{p, ^ 2 ) 7 ^ (0,0) and in Case 3 A 2 7 ^ 0. Due to the presence of arbitrary constants Ai and A 2 , the 
constant e can be assumed to take the values ±1 only. 

The following example shows how to recover the group classification of class ([T]) using the 
results obtained for class m- 

Consider Case 1 of Table 3 extended by the equivalence transformations from i.e., the 
first case presented in Remark 3, were h = Ai(et + k)^, I = X 2 {et + k)^~ . We denote fg{t)dt 
by T. As t = T and 1{T) = {ffg)tfg, we get {f/g)t = X 2 g{t){eT + k) V. Finally, 


f{t) = 


^2g{t) i^T + k) 2 + As^ 

X 2 g{t) (i ln(eT + k)-\- A 3 ) , 


if p / 2, 
if p = 2. 






After re-denoting the constants Xi, i = 1,2,3, it is easy to see that we get Cases 1 and 2 of 
Table 2, respectively. To obtain the corresponding Lie symmetry operators one should make the 

p—2 

change of variables t = T, x = x,u = u + H 2 {£T + k)^~ + pa (resp. u = u + fj ,2 ln(eT + k) + 
for the second case) in the vector fields Xi = dx and X 2 = 2(et -|- K)di + epxdx + e{p — 2)udu- 
It is interesting to note that the images of two distinct inequivalent cases of Lie symmetry 
extensions in class (1) (Cases 1 and 2 of Table 2) belong to the same case of Lie symmetry 
extensions for class (10) (Case 1 of Table 3). 

The other cases are easily treated in the same way. 

5. Conservation laws 

We classify (local) conservation laws of equations from class ([T]) , applying the most direct method 
based on the definition of conservation laws. 

We briefly present necessary notions, specifying them to the relevant case of a single partial 
differential equation C: L = 0 with the two independent variables t, x and the single dependent 
variable u, where L is a differential function, i.e., a smooth function of t, x and derivatives 
of u. A tuple {F, G) of differential functions is a conserved current of L if DiF + DxG = 0 
for any solution of T, where Dt and Dx are the operators of total derivative with respect 
to t and x, respectively. Then the functions F and G are called the density and the flux of 
the conserved current {F,G). Conserved currents {F,G) and {F',G') are equivalent if there 
exist such differential functions F, G and H that F and G vanish for all solutions of C and 
F' = F + F + DxH, G' = G + G — DtH. Elements of the factor-set of conserved currents 
with respect to this equivalence relation are called conservation laws of C. Each conservation 
law F of C admits the representation DtF + DxG = XL for some conserved current {F, G) 
containing in F and a differential function A, which is called a characteristic of the conservation 
law F. Characteristics A and A are equivalent if their difference vanishes on solutions of C. Under 
certain natural condition on T, there is one-to-one correspondence between conservation laws 
and equivalence classes of characteristics m- 

Computing conservation laws for equations from class m, we assume that components of 
conserved currents are of the general form 

F = Fit, X, U, Ut,Ux), G = G{t, X, U, Ut,Ux,Utt, Utx,Uxx)- 

Following [1], it is possible to show that in fact any conservation law of each equation from 
class dU possesses a characteristic of order not greater than two and, therefore, contains 
a conserved current whose density and flux orders are not greater than one and two, respectively. 
The proof of this is quite long and will be a subject of another paper. 

The classihcation of local conservation laws of equations from the class © is as follows. 

Case 0. Each equation from the class admits the “natural” conservation law with the 
constant characteristic A^ = 1. The corresponding density and flux are 

F^=u, G^ = f{t)u + ^g{t)u^ + h{t)utx- 

For general admitted values of the arbitrary elements /, g and h the associated space of 
conservation laws is one-dimensional. 

Case 1. If the arbitrary elements satisfy the equation {{l/h)t/g)t = 0, and, therefore, 

Ht) = J 9{t)<it + p2^ , 

where pi and p 2 are constants with (/ 0 i,/ 02 ) 7 ^ (0,0), then the space of conservation laws of the 
corresponding equation of the form ® is at least two-dimensional. The second basis conservation 


law can be chosen to have the following characteristic, density and flux: 


= S- = /^<*>*) 

^ ~ ~ j (^/(^)^ + ^9{tW + h{t)utx^ , 

where {ljh)tlg = pi = const. 

Using the family of point transformations t = pi f g{t)dt + p 2 , x = x, u = u/pi 
related to the group Gi , we can reduce any equation of this case with pi 7 ^ 0 to the form 
ut + f{t)ux + uux + t~^utxx = 0 (tildes are omitted in the latter equation). 

Case 2. One more case with at least two-dimensional spaces of conservation laws is given by 
the arbitrary elements satisfying the condition {{f/g)t/g)t = 0 , i.e., if 

f{t) = ^<71 jg{t)dt + a 2 ^ g{t), 

where ai and <72 are arbitrary constants. The second basis conservation law can be chosen to 
have the following characteristic, density and flux: 

A3 = - 2aix + 2-utx, - 2aixu - 

g 3 3g[ty 

^ - 2aih{t)xutx - 2aif{t)xu - aig{t)xu^ + 

where W = u + f{t)/g{t), if/g)t/g = ai = const. 

Using the family of point transformations i = ai J g{t)dt + cr 2 , x = aix, u = u related 
to the group GJ", we can reduce any equation of this case with ai 7 ^ 0 to the form 
Ut + tUx + uUx + h{t)utxx = 0 (tildes are omitted in the latter equation). 

Case 3. The maximal dimension of the spaces of conservation laws for equations from class © 
equals three and is reached for the intersection of Cases 1 and 2 , where arbitrary elements satisfy 
the both constraints, {{f / g)t/g)t = 0 and {{l/h)t/g)t = 0. Then for each of the spaces, a basis 
consists of conservation laws with the characteristics A^, A^ and A^ and the conserved currents 
(F^,G^), (F^,G^) and {F^,C^), respectively. The corresponding equation can be reduced to 
the form ut + {(7it + cr 2 )ux + uUx + {pit + P 2 )~^Utxx = 0 by transformation ([5]). The further 
simplification is possible by transformations from the group G^. For example, we can set one of 
the linear combinations ait + a 2 or pit + P 2 to t ii ai ^ d or pi 7 ^ 0 , respectively. 

A well-studied subcase of Case 3 is constituted by constant-coefficient equations, for which 
pi = ai = d [nsiis]. Up to Cl -equivalence any constant-coefficient equation from class © can 
be mapped to the equation ut = uUx + £Uxxt, where s = sgnh = ±1. Then the characteristics 
and the components of the conserved currents of the above basis conservation laws take the form 
(cf. [g p 195]) 


A^ = 1, F^ = u, C^ =- eutx, 

A =eu, = — + ey, - enut^, 

A^ = -|- 2eutx, F^ = G^ = eut — ufx — ^u^utx — 










6. Conclusion 

The aim of the present paper is to enhance and generalize existing results on Lie symmetries and 
conservation laws of variable-coefficient BBM equations of the form ([T|). Comparing the results 
of [7] with those collected in Table 2, we conclude that Lie symmetry extensions presented in [7] 
are particular specifications of Cases 1-5 from Table 2 for certain fixed values of the arbitrary 
element g. For example, there are two cases in the classification list derived in [7] with the 
three-dimensional maximal Lie symmetry algebras (Cases 4 and 10 of Table 1 therein). These 
cases are particular subcases of Case 5 of Table 2 ior g = go = const and g = goe^^. 

The results of Section [5] on conservation laws are quite expectable and, at the same time, are 
not trivial. They naturally generalize well-known results of constant-coefficient BBM equations 
and need the completion of the most significant and tricky part of the proof, which is deriving 
an upper bound for order of conservation laws similarly to 1^. 
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